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Summary 


Wc  present  a  new  formulation  of  the  Taylor  series  for  a  class 
of  functions  which  can  be  expressed  as  the  product  of  a  purely 
angular  and  purely  radial  part,  viz.. 


G(r)  =  YXu(r)f(r) 


in  which  YXjj(r)  is  the  spherical  harmonic  function.  The  formulation 
is  based  upon  the  use  of  the  Fourier  transform  of  the  function  to 
be  expanded.  The  general  expression  for  the  Taylor  series  which  we 
obtain  is  similar  in  form  to  the  Laplace  expansion  for  the  Coulomb 
potential.  Thus,  Taylor  series  can  be  developed  for  arbitrary 
functions  along  lines  which  arc  similar  to  the  multipolar  expansions 
in  the  electromagnetic  theory.  Such  expansions  are  easily  adapted 
to  the  symmetry  of  a  collection  of  sources. 


\ 


l 


1  .  Introduct ion 


Our  purpose  in  this  paper  is  to  present  a  new  and  powerful 
method  of  handling  the  expansion  of  a  class  of  functions  in  a  Taylor 
series.  The  elements  of  this  class  are  simply  those  functions  which 
can  he  expressed  as  a  product  of  a  purely  radial  and  purely  angular 
part :  viz. 9 


CUr)  =  (1.1 

A 

in  which  (r)  is  the  spherical  harmonic  function.  We  proceed 
to  derive  the  general  term  of  the  Taylor  series  to  any  arbitrary 
order.  We  are  able  to  show  that  given  the  Fourier  transform  of 
any  well-behaved  function  (or  generalized  function),  it  is  possible 
to  obtain  the  specific  form  of  an  arbitrary  term  in  the  series. 

These  results  extend  and  generalize  an  initial  effort  which 
was  presented  elsewhere  by  us  (Schmidt,  Pons,  and  McKinley,  1980). 

In  particular,  by  more  laborious  means,  we  derived  specific  formulae 
for  the  first  and  second  order  terms  of  the  Taylor  series.  The 
method  which  we  used  in  that  earlier  work  we  found  not  to  generalize 
easily  to  the  expression  of  a  general  term  for  the  Taylor  series. 

The  results  of  this  paper  contain  exactly  the  results  reported 
previously  (Schmidt,  ct  al . ,  1980). 

We  find  several  reasons  for  undertaking  the  development  of  this 
new  expression  of  the  Taylor  series.  foremost  is  our  need  to 
consider  the  displacement  of  a  particle  about  a  point  of  equilibrium 
which  simultaneously  is  also  a  centre  of  symmetry  for  a  collection 
of  sources  (Schmidt,  et  al.,  1980).  Thus,  we  need  to  consider 
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a  form  of  the  Taylor  series  expansion  of  an  arbitrary  function 
which  exploits  the  symmetries  of  the  system  in  the  same  manner 
as  is  customarily  clone  for  the  electrostatic  potential  through 
the  use  of  the  Laplace  expansion  (£f#,  Jackson ,  1962).  These 
symmetries  are  most  effectively  handled  through  the  use  of  expansions 
which  depend  upon  the  spherical  harmonic  functions.  Hence,  it  is 
natural  to  consider  the  use  of  Fourier  transforms. 

Often,  and  for  molecular  physics  in  particular,  it  is  necessary 
to  consider  approximations  of  complicated  functions.  This  is 
especially  true  of  potential  energy  functions  which  are  derived 
from  molecular  integrals.  The  analysis  of  the  mechanics  of  molecular 
vibrations,  for  example,  commonly  depends  upon  approximations 
which  replace  the  exact  potential  energy  functions  with  harmonic 
replicas.  it  is  universally  known,  of  course,  that  these  harmonic 
approximations  are  adequate  only  when  applied  to  the  analyses 
of  the  molecular  low-lying  vibrational  states. 

It  is  frequently  of  interest  to  examine  higher  order  terms  in 
order  to  assess  the  accuracy  of  the  harmonic  approximation.  In 
some  cases,  higher  order  terms  are  required  in  order  to  establish 
limits  of  stability  for  mechanical  systems.  For  a  relatively  simple, 
complete  function,  it  is  possible  to  extract  cartesian  harmonic  and 
higher  order  terms  as  individual  terms  in  a  Taylor  series.  Realistic 
potential  energy  functions,  however,  often  contain  fairly  complicated 
angular  dependencies.  These  angular  dependencies  reflect  the 
complicated  character  of  the  environment  which  surrounds  a  particular 
particle  (atom  or  molecule)  of  interest  [see,  for  example,  Briels 
(1980)].  For  these  functions,  the  forming  of  the  Taylor  series 
ceases  to  be  a  simple  matter. 
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The  usual  representation  of  the  Taylor  series  presents  the 
determination  of  the  coefficients  as  a  sequence  of  differentiations. 
These  operations  are  applied  to  the  function  the  expansion  of  which 
is  required.  It  is  possible  in  principle  to  make  repeated  use  of 
Rose’s  (1957)  formula  for  the  gradient  to  determine  these  coefficients 
in  terms  of  the  basis  vectors  of  the  spherical  tensors.  In  reality, 
however,  any  actual  attempt  to  apply  that  formula,  even  to  the 
simplest  of  functions,  quickly  becomes  unmanageable.  It  is  tedious 
even  to  determine  the  harmonic  terms  for  the  series. 

In  contrast,  we  are  able  to  obtain  the  coefficients  of  the  Taylor 
series  after  the  evaluation  of  only  a  single  form  of  radial  integral. 
The  strength  of  the  method,  therefore,  is  a  transparent  compactness 
and  flexibility  in  terms  of  ease  of  application.  This  strength  is 
illustrated  by  several  examples  in  a  separate  paper  (McKinley  and 
Schmidt ,198  b) . 

The  outline  of  this  paper  is  the  following.  In  section  2 
we  consider  the  single  centre  expansion  of  a  function  into  a  Taylor 
series.  A  particular  form  of  our  result  replaces  differentiations 
with  integrations.  For  complicated  functions,  which  possess  well- 
defined  Fourier  transforms,  these  integrals  prove  easier  to  consider 
and  evaluate  than  is  the  case  of  the  consideration  of  the  multivariate 
differentiation  of  the  function  directly. 

In  section  3,  based  on  the  use  of  the  results  of  section  2,  we 
present  a  derivation  of  Rose's  (1957)  formula  for  the  gradient 
operation.  A  consistency  between  Rose's  work  and  ours  is  established. 
In  addition,  we  use  the  differential  raising  and  lowering  operators 
for  the  spherical  Bessel  functions  to  re-establish  a  differential 
form  for  the  Taylor  series  coefficients.  This  differentia]  form 
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is  consistent  with  the  angular  dependencies  of  the  integral  form 
and  with  the  symmetries  of  the  spherical  harmonic  functions.  It 
is  clear  from  the  differential  form  obtained  that  in  many  instances 
the  integral  representation  is  indeed  the  simpler  way  to  approach 
the  determination  of  coefficients  in  a  Taylor  series. 

In  section  4  we  consider  the  Taylor  series  as  an  expansion 
about  two  centres.  It  is  possible  to  consider  the  general  develop¬ 
ment  of  any  well-behaved  function  as  an  expansion  about  two  or 
more  centres.  The  Carlson-Rushbrooke  (1950)  expansion  of  the  free- 
space  Green  function  is  such  an  example.  In  section  4  we  establish 
a  parallel  treatment  of  the  Taylor  series.  We  limit  our  considerations, 
however,  only  to  two  centres.  In  this  section  we  also  show  that  as 
the  two  centres  coalesce,  the  resulting  form  is  consistent  with  the 
single  centre  form  of  the  Taylor  series. 

Finally,  in  the  last  section  we  mention  the  limiting  forms 
for  scalar  functions.  Whereas  for  the  general  case,  for  a  general 
angular  dependence,  the  integral  form  of  the  Taylor  series  seems 
to  be  the  simpler  route  to  follow,  for  a  scalar  function,  we  find 
the  resulting  differential  form  is  equally  simple,  if  not  simpler,  to 
use.  Our  form  of  the  Taylor  series  in  its  differential  representation 
for  a  scalar  function,  however,  does  not  hear  a  close  resemblance 
to  the  cartesian  differential  form.  Ou**  form  preserves  the  angular 
dependencies  as  arguments  of  the  Legendre  polynomials.  Differentia¬ 
tion  is  applied  only  to  the  radial  part  of  the  function. 
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2.  The  Taylor  Ser i es  for  the  Expansion  about  one  Cent r e 


We  develop  the  Taylor  series  for  functions  which  are  separable 
into  purely  radial  and  purely  angular  parts,  viz.. 


C(r)  =  F(r)YAu(?) 


(2.1) 


where  (r)  is  the  spherical  harmonic  function.  Still  more  general 
functions  could  be  expressed  as  superpositions  of  such  functions. 

The  function  G(r)  admits  a  Fourier  transform: 


G  (  r ) 


d  3  k  f ( k ) Y  A  ^ (k)exp(-ik*r) 


12.2) 


where  f ( k )  is  given  by 


r  cm 


4  TT  i 


dr  r2F(r)jx(kr) 

0 


(2.3) 


and  j^(x)  is  the  spherical  Bessel  function  of  the  first  kind. 

The  function  at  the  displaced  point  r  +  c  is  given  by  the 
vectorial  Taylor  series  (Band,  1959) 


«:(r  +  c)  =  l  fc-V)n(’,(r).  (2.4) 

n  =  0 


Consider  now  the  general  term  in  this  series  with  G(r)  given  by 
its  F . t  .  : 
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n ! 


tc  -  V) 


(Hr) 


1 


(2tt)  3n ! 


d3k  f t k ) V ^  (k) ( - ic »k)  cxp(-ik*r) 


C~-L_cJ-- fd  3  k  knf(k)Y.  (k)  (c  •  k )  "exp  ( -  i  k  •  r)  . 
C  2 7T )  3  n !  1  Am 


Here,  c • k  is  the  cosine  of  the  angle  between  the  unit  vectors  c 
The  nth  power  of  this  quantity  can  be  expanded  in  Legendre  poly 
nomials  of  the  same  argument  (Morse  and  Feshbaek,  1953) 


A  A  /\  /\ 

(c-k)n  > 


The  coefficients  ^  are  given  by 


A 


2£+l 


nil 


+  J 


dx  x  Pjjfx) 


0 ,  £>n  or  n- £  odd 


(.  2  £  + 1 )  n ! (n- £+1 ) ! !  „  ,  „ 

=  (n-£+l) !  (n  +  g  +  OTT’  *-n  and  n'£  evcn‘ 


We  now*  introduce  the  Rayleigh  expansion, 


exp  ( -  i  k  •  r)  =  \  ( -  i ) l'  ( 2L+ 1 )  j ,  (kr )  p  (k  •  r )  , 


and  the  addition  theorem  for  the  spherical  harmonic  functions, 


f  1^  ( k  •  r )  2I+1  l  Y|,M^k^LMlr 

M 


2.5) 

a nd  k . 


2.6) 


(2.7) 


12.8) 


so  that  we  have 
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exp  ( -  i  k  •  r )  -  \  7i  y  (  -  i  )  'j  j  l  kr  |  Vj  Mlk  IVj  M( 
L  ,  M 


(c-J)"  .  t-  I  Y*  (kivtm(c). 

2 


(2.  11) 


With  these  substitutions,  Equation  (2.5)  becomes 


Jr(c-V)nr.(rl  .  l  (-i)1-*n(2l*n'\lYLM(?)Y,11(0 


L ,  M  ,  m 


x  dk  k2  jdnkknf(k)jL(kr)YXp(k)Y*M(k)Y*m(kl .  (2.12) 


The  angular  integral  is  the  well-known  integral  over  three  spherical 
harmonic  functions  (Rose,  1957) 

j  O-YooUPY u.wi.1  ».1 


2.13) 


where  (L?Mm|Au)  is  the  Clcbsch-Gordan  coefficient.  We  also 

introduce  a  quantity  defined  by  the  remaining  radial  integral: 


InI  (r)  =  — —  dk  kn  +  2f(k) j.  (kr) . 

nL  ( 2tt )  3 ^  o  L 


(2.11) 


Altogether,  we  have  found  the  following  expression  for  the  general 
term  in  the  Taylor  series: 


Jrlc.YI-CUr)  •  KT  I  { -  i ) *‘*nA  (Yj  >j( r ) V  f c  1 

L  ,MJt  ,m 
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x 


ZL+ 1 

(22  +  1) (2X  +  1 ) 


1/2 

(L200 | XO) (L2Mm| Xu) I  , (r) . 


(2.15) 


The  most  frequently  used  cases  will  he  those  of  the  lowest 
order.  For  n=l,  2  can  only  be  1  and  A^  =  l.  We  have  then 


c • VG (r) 


(4u) 3/2c 


l 

L ,  M ,  m 


(•i)L*1YLM(hYlm(ch 


2L  +  1 

n 


Fl) 


i/: 


x  (L100  I  XO)  (LlMml  Ag)IH  (r)  . 


(2.16) 


A  special  application  of  this  is  to  find  an  equilibrium  condition, 
which  can  be  expressed  as  c#VG(r)  =  0,  for  any  c,  or  as 


0=  l  (-i)LYTM(^)(2L+l)1/2(L100|X0)(LlMm|Xu)I1I  (r).  (2.17) 

L,M  iL 

For  n=2,  2  can  be  either  0  or  2  with  A2q=1/3  and  We  have 

then 


3/2 


L,  2  r  /•  s  (4ir)J/i- 
2-(c*V)zG(r)  =  -i - -  c 


*  I  A  (•i)L*2vLMd)V2„(c) 

L  ,M,m 


2L  + 1 


5 (2X  +  1 ) 


x!2L(r) 


(c)I2x(r) 

1/2 

( L  2  0  0 | XO) (L2Mm| Xp) 

(2.18) 


These  lowest  order  terms  are  obviously  easy  to  evaluate,  but  even 
the  general  term  [F.quation  (2.15)]  is  not  intrinsically  more 
difficult.  This  simplicity  traces  back  to  Equation  (2.11)  which  is 
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an  especially  simple  example  of  the  expansion  of  an  invariant  into 
irreducible  forms  (Fano  and  Racah,  1959). 

3.  The  Formula  for  the  Gradient 


In  this  section  we  verify  that  the  first  order  term,  eqn  (2.1b), 
is  equivalent  to  the  usual  gradient  formula  (Rose,  1957). 

We  begin  by  evaluating  the  Cl ebsch-Gordan  coefficient  (1.100|>0) 
in  expression  (2.1b).  This  coefficient  restricts  the  parity  of  L 
as  well  as  limits  its  range.  We  have 


1/2 

(A-1,100|AO)  - 


(A  +  l,  100|  AO) 


A+l]1/2. 


(3.1) 


all  other  values  of  L  yield  zero.  We  now  write 


c*  VG (r)  =  (-i) 


A  (4tt)5/2  ^ 


/T  M,m 

A  +  l) 1/2 


2A+T 


l/: 


(A-l,lMm|Ap)YA_1M(r)Ylm(c) 


2JTJ 


(X  +  MM»|Ap)Yx  +  lMrr)Yln(c)I,fX  +  1(r) 

(5.2) 


We  handle  the  angular  and  radial  factors  separately.  The  standard 
definitions  of  the  vector  spherical  harmonic  functions  are  (rose, 
1957) 


Ta.a  +  i.u^ 


J  H(?)C 

M  ,m 


(3.3) 


where  the  quantities  ^  are  the  basis  vectors  for  the  spherical 
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tensors  (Rose,  1957).  The  scalar  product  of  (7.3)  with  the 
vector  c  i s 


STxXMpfrl  =  J  O^Mm|Xu)Y 

M ,  m 


(3.4) 


=  /GTH*  l  0±llMm|Xu)YA  +  1M(r)Ylm(c) 
M,m  ll 


On  the  other  hand,  for  the  radial  part  we  consider  the 
following.  In  the  definition  of  we  usc  the  standard 

operators  for  lowering  and  raising  the  index  of  the  spherical 
Bessel  functions  (Morse  and  Feshback,  1953,  and  Infold  and  Hull, 
19  51).  Thus, 


(2*)3ln-i(r)  = 


dk  k3f(k).iA_j  (kr) 


dk  k 3  f (k) 


'  0 


T Thfkr)  +  TIT7T  jx(kr) 


|A  +  1  d 
{—  Jr 


dk  k2  f  (k)  j  A  (kr)  . 


(3.5) 


Similarly, 


(2TT)3In  +  1(r)  =  j  dk  k3f(k).ix  +  ]  (kr) 


'X  _ 

d  ] 

r 

d7j  j 

dk  k2  f (k ) j  (kr) . 


3.(0 


At  this  point,  wo  identify  the  inverse  of  eqn  (2.5) 


Hr) 


-  L~lL 


2  7T  2 


dk  k 2  f (k) j . (kr) 


(  3  .  "  ) 
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Wpon  substitution  of  all  of  these  expressions  into  eqn  (3.2), 
we  find  exactly  the  standard  formula  for  the  gradient  operation 
(Rose,  1957) : 


This  new  form  in  terms  of  Fourier  transforms  could  have  been 
anticipated,  because  the  radial  operations  in  eqn  (3.8)  are  uniquely 
those  associated  with  the  spherical  Bessel  functions. 

Fqn  (3.8)  completes  the  connection  with  Rose's  formula  for 
the  gradient.  We  now  extend  the  analysis  given  above  to  the 
consideration  of  an  arbitrary  term  in  the  Taylor  series.  Thus,  we 
consider  the  transformation  of  an  arbitrary  term  from  an  integral 
to  a  differential  form  in  the  radial  component  only. 

For  a  given  function  (l(r),  A  is  specified.  For  a  given  value 
of  the  order  n,  the  values  of  l  also  are  specified.  The  Clebsch- 
flordan  coefficients  fix  the  parity  and  range  of  the  L-values. 

Thus,  it  is  possible  to  see  that 

b  =  A  +  q  (3.9) 

with  the  values  of  q  fixed  and  limited  by  the  conditions  mentioned. 

The  summation  over  in  the  expression  for  an  arbitrary  term  in  the 
Taylor  series  is  rewritten  as 
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(c-VVY,(r)  =  (.411) 


3/2  c 


l  (-n1 


[  (2f.+n  (2e+n  l 


7T7T  Y«.mU  ) 


-  l  l  (-i: 

q  =  0  M 

+  (-i )'q/T 


X+q ) + 1  (X+q £00  |  X  0)  (A+q*Mm  j  Ap)  +qM(r)  In  +q  ( 
r+T(X-q£00 | XO) (X-qHMm|Ay)Yx_qM(r)In>A_q(r)  . 


(3.10) 


In  fact,  the  summations  run  over  q=0,2,4,6...  or  q=l,3,S,... 
depending  upon  tlie  values  of  n  and  X.  The  factor  togetlier  with 

the  Cl ehsch-Cordan  coefficient  automatically  sort  out  the  terms 
which  survive  in  the  summation.  Hence,  there  is  no  need  to  display 
any  specific  form  of  the  sorting  process. 

We  now  use  the  following  lowering  and  raising  operators 


-  c-nq  C*qU-]VJ<k-Y(k,.)) 


J**MCkr)  ■  (‘1)  r  ImrJ  r  lk 


(3. 1  1  ) 


Jx-q<kT> 


rq_''1  (V3v)  q-^CkT"'  } 


(3.121 


to  write 


■  eirr  ^ 


r'x[  dk  kn'q*2f(k).i,(kr)  (3.1J 


(Zit)  31  (r)  -  r<|-X-|U7)V*,|  dk  kn-<'*2f(k)jl(kr)  O. 14) 

f  \  J  J  ft 
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For  a  given  value  of  q  in  1  .  (r)  we  have  adjusted  the 

n  ,  a  7  q 

index  on  the  Ressel  function  to  a  value  of  X.  It  remains  therefore 
only  to  lower  the  magnitude  of  the  exponent  of  K  in  the  integrand 
from  kn  +  "  to  k~.  This  we  do  as  follows.  First,  note  that  from 
the  various  conditions  on  the  parities  of  the  indices,  n-q  will 
always  be  even.  Thus,  we  can  reduce  n-q  to  a  zero  value  in  steps  of 
2.  From  the  eigenvalue  equation  which  the  Ressel  function  i ^ 
satisfies,  we  write 


x2.i  A  (x) 


X  (  A  + 1 ) 


2  d 
37 


U(x) 


(3.15) 


The  right  hand  side  of  this  equation  defines  an  operator  which  lowers 
the  exponent  on  k  (x=kr)  by  2.  The  use  of  this  operator  (n-q)/2 
times  yields  the  desired  result: 


dk  kn~q+2f (k) jx(kr)  = 


=  2  tt  2  i 


X 


X (X  +  l ) 


n-u 


x(x+n 


d  (  r  2  d  X 

37^r  37' 


J  i  v2  3  X 

37u  37' 

n-q 


j  dk  k 2  f  f  k 1 i ^ ( k  r ) 


~T  '  i-n,,, 

r  1  Fir). 


( 3 . lb) 


The  final  expression  for  the  arbitrary  term  in  the  Taylor 
series  is 


n  £ 


±  (cV)nG(r)  =  /4u  l  (-i)n  - 

n-  ~  ~  ~  n>  9,m  [  (2X  +  1  (25  +  1)] 


T77  Yfm(c) 


'll* 

(pn  m 


✓T(TTq7TT(A  +  qff,00|X0)  ( X  +  q PMm  |  X p )  Y ^  +  , ( r )  r 


X^q 


rd 7 


-  A 


/2  ( X -c,)  + 1  ( > -q  POO  |  AO)  fA-qtMm|Ap)Yx.qM(r)rq'X_1[7j?]  V  +  1 


1 
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*  A(A+1) 


ra_nF(r) 


f  3. 17) 


ision  o!  a  Function  in  a  Tavlor  Series  about  two  Centres 


lor  some  applications  it  is  useful  to  consider  the  function  of 
ei|i)  (2.1 )  to  depend  upon  two  points  n  and  r2  ,  either  of  which 
may  he  displaced  independently.  We  set 

r  =  r,  -  r 2  (1.1) 

and  redefine 


Cfr,  ,  r 2  )  =  Y^l  r)  F  ( r)  .  (1  .  .1 ) 

I’ll  is  function  still  possesses  a  Fourier  transform  as  expressed  by 
eqn  (2.2),  (2.3)  and  (3.7).  Before  we  evaluate  this  function  at 
displaced  points,  however,  we  present  an  alternative  Fourier  transform 
which  emphasizes  the  two-centre  character. 

Into  the  expression  for  the  Fourier  transform 


(•  (r  i  ,  r 2 )  =  (2n) 


d 3  k  f(k)YA^(k)exp(-ik*ri)cxp(ik*r2) 


(4.3) 


we  substitute  Luw  Rayleigh  expansions  of  the  form  (2.10).  The 
angular  integrations  over  the  three  spherical  harmonic  functions 
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G  t  r  i ,  r  2 ) 


m  O 

'I. 


dkk 


,t  2 


dBkf  (k)YXll(k)  & 

4 7T  L i  ,  M i  ,  L2  ,  M2 


XjL1^rOVLiMi(;i)Y*iMi(k)jL2(kr2)Y*2M2(;2)VL2M2(k) 


1  E  iL2_LlYL|Ml^^YL2M2^) 


tr3/2  L1M1L2M2 


f  CZL,1^!)  C2L2^l)j  1/2(Lil2oo|  AO)  (Li  L2MjM2  |  Ay) 


r 

xj  dk  k2 f (k) jLi (kri ) jL2 (kr2) 


(4.4) 


In  order  for  eqn  (4.4)  to  be  consistent  with  eqn  (3.7)  and  (4.2) 
for  arbitrary  functions,  we  require  the  following  double  addition 
formula  for  spherical  harmonic  and  Bessel  functions: 


YAMtr)jA(k1-)  =  (4ir) 
f(2Li+l) (2L2+1)1 

- zrn 


1/2  £  iA  +  L2 -L ! (LlLa00| A0) (L1L2MiM2 I  An) 

LiMjLzMz 

1/2 

Yl.1M1(ri):lLx(kri)YL2M2(r2)jL2Clcr2)  •  (4-5) 


A  similar  form  for  real  spherical  harmonic  functions  has  been 
derived  by  Johnson  (1973),  It  is  possible  to  show  easily  that  the 
right  hand  side  of  (4.5)  reduces  to  the  form  of  the  left  hand  side 
when  r2  =  0  and  r=r i  .  Additionally,  when  A  =  y  =  0,  the  right  hand  side 

A 

of  (4.5)  can  be  shown  to  reduce  to  Y0o(r)jo(kr)  with  the  use  of  the 
appropriate  Clcbsch-Gordan  coefficients  and  the  addition  theorem 
for  the  spherical  harmonic  functions.  Thus,  the  two-centre  expansion 
is  consistent  with  the  more  familiar  single-centre  forms. 

We  proceed  to  evaluate  the  two -cent re  function  at  the  displaced 
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points  rj+a  and  r2+b.  We  have  as  the  extended  vectorial  Taylor 
series,  the  expansion 

+  b • V2 ) nG (r  j , r  2 ) .  (4.6) 

is  to  be  evaluated  with  G(ri,r2) 
expressed  by  its  Fourier  transform  [eqn  (4.3)].  We  have 

+  b-y2)nG(r,  ,r2)  =  (2tt)  "5Jd3k  f(k)YAji(k) 

x  -pj-(-ia*k  +  ib*k]n  exp  ( -  ik  •  r  i )  exp  (i  k  •  r  2 )  .  (4.7) 

In  this  expression  we  use  the  binomial  expansion  and  the  expansion 
in  irreducible  forms  (2.11)  to  replace  the  quantity  in  the  brackets: 


oo 

G  ( r  i  +  a  ,  r  2  +  b )  =  £  -i-(a-V! 

~  ~  ^  "w  _ n  1 1  .  ~  ~ 


n=0 


The  general  term  in  this  series 


[a-k  -  b - k ] n  =  £(") (a-k)n*q(-b-k)q 

q  1 


y  fn>.n  n-qr  (4tt)2 

0  \  VK  d  (  bJ  XHT+Ij  ( 2 £ 2  +  1 )” 

q , £ ,  ,  £  2  m i , m  2  '  1 


xAn-q,tlAq,«.2YK.imi  ^a^h1mdk)VP.2m2  ?2m,  ’ 


(4.8) 


The  coupling  rule  for  the  spherical  harmonic  functions  is  now  used 
to  reduce  the  number  of  these  functions  in  the  product  (Rose,  1 9 S 7 1  : 


1/2  U.  ooo  I  »0) 


*Ui  e2mijn2  I  ^n0  Y  £m  (  k )  . 


(  4 . 9  ) 
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The  substitution  of  two  Rayleigh  expansions,  and  the  use  again  of  the 
coupling  rule  together  with  eqn  (4.8)  and  (4.9)  leads  to 


~r[a-V,  +  b-V2]nG(r,,r2)  =  I 

q,£i,£2,mi,m2>Li,L2,Mi,M2,L,M 


x(_i)n+L1-L2(n)an-q(-b)q 


An-qf£,Aqf£,YtimiWYt,»,^YLiM^riJ 


xYL2M2  l'r2')  TTT7+T}  fziT+li  (2£+TTT2L+1)  ,  ( «. i £.2 00 1  £0)  (fci£2niim2 1  Cm) 


x(LjL200|L0)  (L1L2M1M2 |LM)J  dk  kn+2 j L ^  (kr 1 ) j j ^ (kr2 ) 


x.r47idVxw^)<M^)v;m(k). 


(4.10) 


The  angular  integral  is  evaluated  as  before  (cf.,  eqn  (2.13)),  and 
we  introduce  another  quantity  defined  by  the  remaining  radial  integral : 


f 00  ~ 

J,  ,  (r  1  ,  r2)  =  dk  kn+2f(k)j,  (kn)j,  (kr2) 

nLil,2  -  ~  ( 2 tt )  3  J0  Li  L2 


(4.11) 


Thus,  we  obtain  the  following  expression  for  the  general  term  in  the 
Taylor  expansion  of  a  two-centre  function: 


Hit?'!1  +  b*V2]nG(ri  ,r2)  =  -  I 

n !  q  ,  ^  1  ,  ^  2  ,  ni  1  ,  m  2  ,  L  j  ,  L  2  ,  M  j  ,  M  2  ,  L  ,  M 

n-i)n*L‘-L!(")0n-'>eb)l'An-q,slAq,tP8l„1O)'-^mj  (lhv,lMl  in) 


C1'^00!10’  Cl.l.ooIlo) 


X(U00|  AO)  (£i£2Hin2|£m)  (LiL2M|M2  |  LM)  (JtLmM  |  Ap).JnJ  ,  (r2 ,r2)  .  (4.12 
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Alternate  expressions  lor  the  expansion  about  two  centres  are 
possible.  We  present  the  formula  for  the  general  term  of  the  Taylor 
series  for  one  such  modification.  In  particular,  when  the  vectorial 
difference  ri  -  i'2  =  R,  a  vector  drawn  between  two  points  of  reference, 
is  a  constant  quantity,  then  the  expansion  can  be  considered  in  terms 
of  displacements  about  the  two  ends  of  the  vector  R.  The  formula 
for  the  general  term  is 


i,[a*Vi  +  b-VzAlip  = 


l 

q,2. 1  a 2  ,m.  ,m 2  ,  2,,I.,M 


II  ~  CJ  ,  i  q  ,  2.  2 


21. +  1 

1(2 2. i+1)  (2t2  +  lTr 2 A'+Tj' 


1/. 


(2 , 2  20()  I  20)(  21.00 1  AO)  ( ti«.2mim2  |  £m)  ( 2  l.mM  |  A(i)  Y  (;i)Y.  (b) 

V  i  III  i  «  2  m  X 


(1.131 


in  which  I  ^  (R)  is  defined  by  cqn  (2,14).  The  analysis  which  leads 
to  eijn  (4.15)  parallels  the  proceeding  analyses. 


5.  Discussion 


Many  specific  examples  of  functions  which  can  he  expanded 
usefully  in  Taylor  scries  easily  spring  to  mind.  Several  examples 
which  can  be  derived  from  the  Yukawa  potential  arc  considered  in  a 
companion  paper  (McKinley  and  Schmidt,  198  a). 

Wc  conclude  this  paper  by  showing  the  limiting  forms  for  the 
Taylor  scries  when  the  function  to  be  expanded  is  scalar. 


1 


The  scalar  function  (i(r)  is  expressed  as 


C(r)  =  Y00(r) [/4?F(r)) 


=  I;(r). 


The  Fourier  transform  of  this  function  is  Y00(k)f(.k)  with  f(k) 
given  by 


f  (.  k )  =  dr  r  2  F ( r ) j  0 (kr ) . 


(5.1) 


(5.2) 


In  terms  of  this  form  of  scalar  function,  the  general  term  for  the 
Tavlor  series  is 


ff!  l£-?)n«Cr>  -  (4i.)3/2(c"/n!)  I  C-0"*\pta(OVLM(r> 

i ,  m ,  L ,  M 

k / ( 1\.*XTC2J*TJ ( L ft 0 0  |  0 0 )  (L£Mm  |  0 0 )  1  ,  (r) 

,  ,  .3/2,  n ,  n  v  ,  ^n+£.  0-0  YS,-m^r^  dm^  .  ,  , 

*  (4n)  (.c  /n!)  \  (-j)  A  # - lnP(r) 


=  /4i(cn/n! )^(-i)n+£An£P£(r.c) ln£(r) 


(5.3) 


This  quantity  is  evidently  a  generalization  of  the  expansion  of  the 
Coulomb  potential  in  a  Legendre  polnomial  series: 


r +sf 1  =  ‘Vr*c)- 


The  general  differential  form  also  assumes  a  simpler  appearance 
in  the  scalar  limit.  We  need  only  consider  the  I n ^ f r )  quantity  in 
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eqn  (5.3).  I  rom  the  definition  of  1  ^(r)  wo  have 


dk  kn+2f(kj jf fkr : 


=  1  If 


5/2(_1)(n+Jl)/2r£ 


d 

r  Jr 


v. (d/dr ) "  Tl’(  r ) 


(5.5) 


following  Todd ,  Kay  and  Silvcrstono  (1970),  wo  can  show  that 


d 

rTr 


r  1  ( d  /  d  r ) n  *  *■  r  F  ( r ) 


y  !  -q 

q  =  0  (2-qT!  (2q)  !  ! 


n  -  q  .  d 
T  dr 


(d/dr)n'q'1l:(r) 


(5.0 .) 


Thus  , 


I  A  r  ) 
n  v1  v 


i  nm^£)/2  f  CiiiluiH)  i 

~  ‘  y  1  >  L  To  -ni  i  r  ->,n  i 

^  4  7i  q  =  0 


fFtriTnirnT  r 


n  -  q  d 

r  J7 


(d/dr)n_q"lF(r). 


(5.7) 


lor  example,  the  second  order  term  in  the  Taylor  series  assumes 
a  s  i  in  p  1  e  f  o  nil :  vie., 


55 


1 


Zl 


2  dl  ,  d2f 
r  J:  dr’ 


ZV  :  ( r • c ) 
.*>  * 


1  dl 
r  d  r 


d/1- 
d  r?J 


(  5 . 8  ) 


This  form  is  particularly  useful  for  a  number  of  applications  which 


involve  scalar  i’unct  ions. 
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